Lack of self-average in weakly disordered one 
dimensional systems 



A. Crisantif, G. Paladin^, M. Serva§ and A. Vulpianif 

t Dipartimento di Fisica, Universita 'La Sapienza' 
P.le Moro 5, 1-00185, Roma, Italy 

I Dipartimento di Fisica, Universita dell'Aquila 
1-67100 Coppito, L'Aquila, Italy 

§ Dipartimento di Matematica, Universita dell'Aquila 

1-67100 Coppito, L'Aquila, Italy 

PACS number(s): 75.10.Nr, 05.50.-fq, 02.50.-fs 



1 



Abstract. We introduce a one dimensional disordered Ising model which at zero 
temperature is characterized by a non-trivial, non-self-averaging, overlap probability 
distribution when the impurity concentration vanishes in the thermodynamic limit. 
The form of the distribution can be calculated analytically for any realization of 
disorder. For non-zero impurity concentration the distribution becomes a self- 
averaging delta function centered on a value which can be estimated by the product 
of appropriate transfer matrices. 
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1. Introduction 



Disordered systems attracted much work in the last years. One of the main features 
of these systems is the large number of locally stable states. As a consequence, a 
natural characterization of the equilibrium states is through the distribution of their 
mutual distance, that is through the probability distribution P{q) of the overlap q 
between pairs of states [1]. In general, the organization of the equilibrium states 
depends on the realization of disorder, so that P{q) depends on the disorder, even 
after the thermodynamic limit has been taken. For instance, this lost of self-averaging 
is observed in the Sherrington-Kirkpatrick model of spin glasses [2] where it can be 
related to a symmetry breaking in the replica space [1] . A big effort has been devoted 
to the search of simpler systems which share the main features of the Sherrington- 
Kirkpatrick model. The most celebrated among them is the random energy model 
[3] now used in very different contexts such as biological evolution, neural networks, 
polymers, and so on. On the other hand, it is commonly assumed that one dimensional 
systems lack for such a rich behaviour. In this paper we introduce a disordered Ising 
model which, in spite of being a one dimensional model with short range interactions, 
exhibits at zero temperature a non-trivial, non-self-averaging, P{q) if the impurity 
concentration vanishes in the thermodynamic limit. This latter condition is crucial, 
since if the disorder concentration does not vanish in the thermodynamic limit, P{q) 
becomes a self-averaging delta function, even at zero temperature. This is the main 
result of the paper which allows one to get a deeper understanding of the mechanisms 
underlying the appearance of non-trivial ergodicity breaking phenomena in disordered 
systems. 

The paper is organized as follows. In Sect. 2, we introduce the model and the 
thermodynamic quantities of interest. In Sect. 3, we study the model in the weakly 
disorder limit, where it is possible to calculate the form of the overlap distribution P{q) 
for any disorder realization. In Sect. 4, the model is analyzed for non-zero impurity 
concentration where P{q) is shown to be a self-averaging delta function. The value of 
the overlap where the delta function is centered can be obtained by means of random 
transfer matrices. 

In Appendix 1, we discuss an extension of the Landau theorem on the absence of 
phase transitions in one dimension to the overlap probability distribution. In Appendix 
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2, we give a brief description of the numerical algorithm used in Sect. 4 to compute 
thermodynamical quantities via product of transfer matrices [4] . 

2. The system 

One of the main reasons for the existence of a large number of equilibrium states in 
statistical systems is the presence of "frustration". In general, a disordered system 
cannot satisfy all the constraints originated from the competing effects of disorder and 
interactions. We introduce a one dimensional model which captures the main features 
of disordered systems but, at the same time, is simple enough for an analytical study. 
The model is a random field Ising chain of N spins described by the Hamiltonian 

N 

H = -^[aiai+i + hiai] (2.1) 

with 

hi^Vi- Vi+i (2-2) 

where the rji are random independent variables which take the value +1 and —1 with 
probability p and (1 — p), respectively. As a consequence, the magnetic field on a site 
assumes the value hi = ±2 with probability p{l — p), and hi = with probability 
+ (1 — p)^. By changing the value of p, we can modulate the disorder in the model. 
The pure ferromagnetic system is obtained for p = 1 or 0, yielding hi = 0\fi. We can 
also consider a pure system with a finite number / of impurities, so that the impurity 
concentration vanishes in the thermodynamic limit, weakly disordered system. The 
integer number / is a random variable distributed according to P{I) = exp(— c)//! 
where c = P{I) I is its mean value. The weakly disorder limit is achieved by 
taking rji = 1 with probability p = 1 — c/N, and hence p — > 1 when the number of 
spins N ^ oo. 

For any value of p the local magnetic field has a certain degree of spatial correlation, 
all field realizations are made by a sequence of strings 

2, O-- - 0,-2, O-- - 02 

where the number of zeros between two successive 2 and —2 or —2 and 2 is a random 
non-negative variable. 
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The overlap between two spin configurations 'a' and '/?' for the same disorder 
realization is defined as 

1 ^ 

This can be usefully computed by introducing two identical replicas of the system. 
The equilibrium probability distribution of qa/s is, therefore, given by 

Piq; N) = Z^^Y.Y1 '^(^ - ^12) n ^^P(-^ (2.3) 

ctI 0-2 a=l,2 

where and H'^ are the Hamiltonians (2.1) of the two replicas and /? = 1/T. The 
normalization factor is the single replica partition function, which can be defined 
in terms of the product of random transfer matrices as 

N N 

= XI n /^[^i^i+i + - ^»)] = Tr n (2.4) 
iji=±ii=i i=i 

where, for this model. 



Aj = with probability p 



e-3/3 



g/3 g-3/3' 



(2.5) 



Aj = I I with probability 1 — p- 



In general P{q; N) depends on the particular sequence of random fields, even in the 
thermodynamic limit N ^ 00. However, one can show that in one dimensional systems 
with short range interactions for any finite temperature P{q; N) is self-averaging, and 
hence P(g; N) converges towards a well defined function P{q) as N ^ 00. This is no 
more true at zero temperature, where the Landau theorem on the absence of a phase 
transitions does not hold, and one could have an ergodicity breaking [see Appendix 
1]. The zero temperature limit of this model, with p — 1/2, was studied numerically 
in Ref. 5 by extracting the diverging part in (2.5), i.e. by writing A^ = Rj with 





= k . or 11 (2.6) 



and using the matrices in the product (2.4). The non-zero value of the zero 
temperature entropy of the model found in Ref. 5 indicates that frustration plays 
an important role. 
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3. Ergodicity breaking in the weakly disordered system 

The weakly disordered system has a great theoretical relevance, because it exhibits a 
non-trivial ergodicity breaking at zero temperature. To our knowledge, this is the first 
disordered one dimensional model where the overlap probability distribution P{q) is 
a non-self-averaging smooth function of q. 

It has to be noted that the pure system, with hi — Vi, also exhibits an ergodicity 
breaking. The overlap between two equilibrium configurations should be zero for any 
finite temperature, so that P{q) = S{q). On the other hand, at zero temperature the 
overlap is = ±1 and P{q) = [5{q — 1) + S{q + l)]/2. This is a quite trivial ergodicity 
breaking due to the zero temperature phase transition in the Ising model with zero 
field. It simply reflects the breaking of the "up-down" symmetry of the model. 

The weakly disordered systems is defined by assuming r]i = 1 with probability 
p — 1 — c/N, where N is the number of spins. As a consequence, there is a finite 
average number c of impurities in the thermodynamic limit, though p = 1. 

Any field realization of the weakly disordered system obeys some simple rules. A 
sequence {rjk-i = 1, Vk = ~1) Vk+i = 1} h^is probability p^ (1 — p) ~ c/N while a 
sequence {r/fe-i = 1, r]k = —1, Vk+i = —1} has probability — p)^ ~ (c/N)'^, and 
the latter event can be neglected in the limit of large N. It follows that if hi = 2, then 
= —2. Therefore in the thermodynamic limit, a field configuration is given, with 
probability one, by "islands" of zero hi of arbitrary length separated by / interfaces 
made of the pair of random fields hj = 2 and /ij+i = —2. 

Within this picture, we can determine the overlap probability distribution by 
combinatorial considerations. Let us consider the case of / = 1 impurity, that is 
particularly simple since there is only one field interface. Assuming periodic boundary 
conditions, the interface can be moved to the boundary, so we can take hi = —2, 
hN = '2- and hi = for i = 2, . . . N — 1. The system is thus equivalent to a pure 
Ising model without external field and with fixed boundary conditions cxi = —1 and 
(Tat = +1. At zero temperature, such a system is frustrated because the spins should 
be alligned with ai on one side, and with ajy on the other one. Consequently, there 
is a spin fiip in an arbitrary point 1 < L < N. The ground state is then made of spin 
configurations with the first L spins down and the remaining N — L spins up. The 
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overlap between two of such spin configurations with spin-flip points L and L' is 

g(L,L') = l-2^^^. (3.1) 

The probabihty distribution is obtained by counting how many times the overlap is 
equal to g, i.e. 

^ N N 

^i(^) = ]^E E^(^-^(^'^')) (3.2) 

L=l L' = l 

since all values of L and V are cquiprobablc. In the limit of large the sum can be 
estimated, by defining the new variables £ — L/N and x — \L — L'\/N, as 

[ d£ [ dx 6{q -l + 2x). (3.3) 
JO Jo 

The integration can be easily performed and yields 

Pi(g) - de e{q-l + 2£) = ^ (3.4) 

It is worth stressing that Pi{q) is the same for all disorder configurations, due to the 
cyclic property of the trace, though it is not a delta function. The average overlap 
q = f Pi{q) q dq = 1 /3 is different from the most probable value q^ip = 1 

Consider now the case of / = 2 impurities, which is qualitatively different from 
/ = 1. Each disorder realization has two islands of zero fields. Assuming periodic 
boundary conditions, one interface can still be moved to the boundaries, i.e. hi = —2, 
Hn = 2. The other must be placed on an arbitrary site j = xN, i.e. hj = 2, 
hj+i = —2. In the limit of large N, the systems with 7=1 and 1 = 2 have the same 
energy but different entropic factors since for / = 2 there are N allowed positions j 
for the second interface. A ground state configuration is individuated by the position 
of the interface, i.e. by the particular value of x in the interval [0, 1]. 

We can regard the system with 1 = 2 impurities as the superposition of two 
systems with 1 = 1 impurity, made of x N and of {1 — x) N spins, respectively. The 
overlap between two ground states of the same disorder realization x is then given by 
the weighted sum of the overlaps of the two systems with 1=1, 

q{x) = xqi + {1 - x) q2 (3.5) 

where the qi are distributed according to 

Pl{<l^) = ^ (3.6) 



At difference with the case / = 1, for / = 2 the overlap depends on the particular 
disorder realization, implying that the overlap probability distribution ^2(9) is not 
self-averaging. For any value of x, ^2(9) is given by 

P2{q;x) = j dq^ j dq2PMPM5{q-q{x)). (3.7) 

By performing the first integral, we get 

P2{.q] = \ j dg2 [{x + q) + q2{q-l + 2x) - ql{l-x)] (3.8) 



where 



q — x\ ( q ~\~ X 

a = max ( —1 , , b = min 1 , ) (3.9) 

1 — X J \ 1 — X 



The overlap probability (3.8) is shown in Fig. 1 for different positions of the second 
impurity (x = 0.1, x = 0.2 and x = 0.5). For comparison, we also report the average 
overlap probability distribution ^2(9)- 

For I > 2 impurities, it is immediate to repeat the above argument. Thus, the 
overlap probability distribution of the system with / field interfaces located on the 
sites xi N, . . . , xi-i N, N is given by 

Pi{q;xi,X2,...,Xi-i) ^ dqi--- dqj 6 I q Xjqj \ Y[Pi{qj) (3.10) 

J-^ -^-1 y j=i J j=i 

with xi = 1 — Yli~^ ^j- This expression takes a simple form in the Fourier space. Let 
us define the Fourier transform of the overlap probability distribution 

/•oo 

P{u;) = / e'^'iP{q)dq (3.11) 



where P{q) = for |g| > 1. In the case / = 1, one can easily see that 

(1 + g) 



Pi 



= (3..) 



Inserting into (3.10) the integral representation of the delta function, we have for / > 1 
impurities we have 



/OO 1 pi /»1 ^ 

— e-'-^ j dq^--- j ^dqi\[ P,{qi) e'" 

^e--^n^i('^^^)- (3.13) 



i=l 
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Therefore, the Fourier transform of the overlap probabihty distribution is given by the 
convolution 



Pi{<jj;xi,X2,...,xi-i)^YlP^i'^^i) with^Xi = l (3.14) 

i=i i=i 

where 

sin(a;a;j) . cos(a;a;^) . sm{uj Xi) 

Pi{u;xi) = 1 h 1 rir- (3.15) 

U> Xi UJXi (WXj)^ 

For any number / > 1 of impurities, the overlap probability distribution depends 
on the particular disorder configuration individuated by the position of the / — 1 field 
interfaces, i.e by the sequence {xi, X2, ■ ■ • xj-i}. The disorder average is, therefore, 
obtained by averaging over the probability of the sequence V{{xi,X2, ■ ■ ■ xj-i}). In our 
model the Xi are independent random variables uniformly distributed in the interval 
[0, 1], so that 

7-1 

V{{xi, X2,--- xi_i}) = H V{xi) (3.16) 

i=l 

where V{x) = 1 for a; e [0, 1] and V{x) — otherwise. 

Since the Xi are independent, in the limit of large / the average overlap probability 
distribution tends to the most probable value Pc{oj) obtained for Xi = l/I, 



P,{io) = [P.iio)]" (3.17) 

since for large / the fluctuations about the mean value c are negligible. For the 
central limit theorem, P^ is a gaussian with mean value q — 1/3 and variance 



(q — 1/3)^ = 2/(9c). The gaussian form very quickly becomes a good approximation, 
as shown in Fig. 2 where Psiq) (the most probable overlap probability with impurities 
at xi = 1/3 and a;2 = 2/3) as well as the average P3{q) are very close to to the gaussian 
with mean value 1/3 and variance 2/27. 

The limit c — > oo, but c/N — > 0, corresponds to vanishing impurity concentration 
withp ^ 1. As a consequence, at zero temperature the overlap probability distribution 
isP(g) =d{q-l/3) for p ^ 1, while it is a double delta P(q) = [5{q-l)+5{q + l)]/2 
for p = 1, i.e. the pure system. 

The weakly disordered system exhibits a non-trivial new phase between these two 
extreme behaviours. For finite temperature, the overlap vanishes both in the pure and 
in the weakly disordered system leading to P{q) = S{q). 
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4. The disordered system 



In this section we discuss the system with non-zero impurity concentration, i.e. with 
hmjv^oo c/N > and p 1 or 0. Unhke the weakly disordered system, in this case we 
cannot perform simple analytical calculations. However, we can employ the transfer 
matrix approach for a numerical calculation [4] . 
The study of (2.3) is equivalent to that of [5,6] 

•A/'ivH^EE^'^"'" n ^^^''^ (4-1) 

ctI ct2 a=l,2 

which can be regarded as the partition function of a system made of two interacting 
replicas, with Hamiltonian + + uNqi2. The term ujN is the macroscopic 
coupling between the two replicas. Only if w = 0, or if the overlap qi2 = 0, A/'7v(<^) 
becomes the square of the usual partition function (2.4). 

The advantage of (4.1) stems from the possibility to write it as a product of random 
transfer matrices [4,5]. The maximum Lyapunov exponent r(a;) of the product yields 
the average of the logarithm of j\/jv(<^), 

N'Niu;) oc exp[Arr(a;)] for iV > 1. (4.2) 

There is no average over the disorder in the l.h.s. of (4.2) since the Oseledec theorem 
[7] ensures that the Lyapunov exponent is a non-random quantity for A?" —> 00, i.e. 
it has the same value for almost all realizations of disorder, a part a set of zero 
probability measure. This is not in contradiction with the fact that the overlap 
probability distribution could be not self-averaging. Indeed, it is possible to prove 
[5] that the self-averaging of the Lyapunov exponent only implies that the extrema 
of the support of P{q) are the same for all configurations of disorder. Consider the 
quantity 

rN{uj) = lnJ e^'^«P(g;iV) dg (4.3) 

where P{q; N) is the overlap probability distribution of a system of size N, and in 
general depends on the realization of disorder. The Lyapunov exponent is related to 
the thermodynamic limit of 

^r^,(c.) = lln y'\^-^P(g;iV)dg (4.4) 
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In the limit of large N, we can insert into (4.4) the saddle point estimate obtaining 
for small uj 



In general, r(a;) is a non-linear function of uj. 

A reasonable ansatz on the finite N corrections to the asymptotic form of P{q) for 
a given disorder realization is 



In this case, the saddle point estimate of (4.4) gives the Lyapunov exponent as the 
Legendre transform of the convex envelope s{q) of S{q): 



As consequence of the Oseledec theorem, (4.6) shows that the envelope s{q) is self- 
averaging in systems with short range interactions. By definition s{q) = for 
q e [9min, Q'max] and s{q) > for q-max < 9 < 1 and -1 < q < qmin- Moreover, 
for large |a;| the saddle point is given by = ±1 so that the asymptotic behavior is 
r(a;) ~ C(-i-) ± oj where C(-i-) are constants. 

This result could appear of rather limited utility, but it assumes a great 
importance when considered as a mark of a replica symmetry breaking. Indeed, if 
limAT^oo P{q, N) = S{q — q), the derivative dr(a;)/du; at = does exist and is equal 
to q. This is the case in one-dimensional systems when T 0. On the other hand, if 
P{q) differs from a delta function it implies a non-differentiable T{(jj) at a; = 0. 

No information is lost if P{q; oo) is a delta function, implying q^i^ = qmax = Q, 
and s{q) has only one zero at q = q. 

The theoretical relevance of this result follows from the possibility of estimating 
the Lyapunov exponent either by a direct numerical calculation or by many analytic 
methods such as weak disorder expansions [8], cycle expansions of appropriate zeta 
functions [9], microcanonical tricks [10-11], and so on. 

This approach has been used in Ref. 5 to evaluate the average overlap q for 
the model with p = 1/2 for both zero and finite temperature. In this reference q 
was obtained from the numerical derivative of r(a;) at a; = 0. Moreover, the zero 
temperature case was done by using an ad hoc trick to extract the diverging part. 




(4.5) 



P{q; N) ~ P{q; oo) + A e'^^'^^ 



r(a;) = max [quj — s{q)] 



(4.6) 
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Here we report some results on the model with < p < 1/2 for both zero and finite 
temperature obtained by using a different method [4] which avoids numerical estimates 
of derivatives. As we discuss in Appendix 2, the entropy s{q) defined in (4.6) is the 



Legendre transform of (l/N)hiJ^N{u). The quantity ^/^{lv) can be written as the 
weighted sum of exp{Nu>q) over the states with a given value of q, 

Mv(a;) = 5^Mv(g)e^"'^. (4.7) 

g 

For large N the saddle point estimate yields 

1 1 

-s{q) = — In J^N{q) = j;^lri J^N{u;)-u;q (4.8) 

where q = q{Lij) is the overlap selected by the chosen value of u. All the quantities on 
the r.h.s. of (4.8) can be computed directly by means of products of suitable transfer 
matrices. A brief description of the method is in Appendix 1. We report the numerical 
results for the average overlap as function of the temperature T — at p = 1/2 in 
Fig. 3 and as function of p at different temperatures in Fig. 4. We have also apply 
our method to the direct computation of the second derivatives, as one can see in Fig. 
5 which gives the behaviour of the spin glass susceptibility 



xsG = {{q-q) ) 

as function of the temperature at p = 1/2. 



u)=0 



5. Conclusions 

We have shown that in one dimensional systems with short range disorder, ergodicity 
breaking cannot appear at finite temperature. In generic situations the overlap 
between two equilibrium configurations does not vanish and has a self-averaging value 
which can be computed by products of random transfer matrices. 

To describe the zero temperature phase transition, a new one dimensional model 
has been introduced. It is complex enough to be a good laboratory for testing ideas 
which can be useful in more realistic situations. In the weak disorder regime, our model 
exhibits an overlap probability P{q) which is a smooth function of q. It does not self- 
averages in the thermodynamic limit but depends on the disorder configuration. The 
P{q) has an analytic form which can be computed in a transparent way without using 
replica trick or other indirect methods. It is an open problem to find analogous models 
in two dimensions. 
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Appendix 1. 

In this appendix, we show that the overlap probabihty distribution of a one 
dimensional random system with short range interactions is a self-averaging delta 
function at any temperature different from zero. 

Consider a one dimensional random Ising model with short range interactions. 
Because of the well known Landau theorem on the absence of phase transitions in 
one dimension, in any realization of the random couplings and/or fields, the spin-spin 
correlation should have the exponential decay 

{(Ti (Ti+r) - (cTi) (cTi+r) ~ e"''/^ (^11) 

where (•) denotes the thermal average with the Gibbs measure exp(— H) for a given 
disorder realization. Since in a typical realization of disorder is the difference of the 
Lyapunov exponent of the product of random transfer matrices related to the model, 
(Al) also follows from the Perron theorem assuring that there is no degeneracy of the 
eigenvalues in the spectrum of matrices with positive non-zero elements. The value of 
^ can depend on the sites i, i + r and on the disorder realization. Nevertheless, (Al) 
implies that the are (exponentially) independent random variables with respect to 
the Gibbs measure. 

The overlap between two spin configurations and cr^ is 

As the ai are independent, we can apply the law of large numbers so that in the 
thermodynamic limit with probability one all spin configurations related to a disorder 
realization have the same overlap, i.e. limAr^oo —>■ Q and the overlap probability 
distribution P{q) converges towards a delta function. 

We should still prove that q has the same value for all disorder realizations. This 
can be seen by means of a very general result: the support of the overlap probability 
distribution P{q) is self-averaging. Consequently, if P{q) is a delta function for a 
disorder realization, it should be a delta function centered on the same value q, for 
all the other realizations. In conclusion, for one dimensional systems with short range 
interactions one has in the thermodynamic limit for any non-zero temperature 

P{q)^5{q-q). 

for almost all disorder realizations, a part a set of zero probability measure. 
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Appendix 2. 



In this appendix we briefly review the numerical method used in Sect. 4 to evaluate 
the average overlap probability distribution. We only describe the finite temperature 
method. We shall denote by Si = (aj^af) the pair of spin at the same site in the 
two replicas. The average of the logarithm of A/Ar(^)j equation (4.1), can be obtained 
from the vector iVj(/S'j_|_i) which satisfy the recursion relation 

Ar,(5,+i) = ^exp(-/3ej - /3e^,+ujala^) N'^_AS^) (B.l) 

where ef = — [cr^^i + hi\af is the single spin energy in the site i of the replica a, and 

J2 KiSi+i) = 1. {B.2) 

From (4.1) follows the initial condition 

By denoting with rii the sum of the elements of iVj(5j+i), the logarithm of A/jv(<^) 
averaged over the disorder is 

I I ^ 

i.e. the Lyapunov exponent of the product. 

The value of q can be obtained directly from the derivative of the vector Ni{Si^i) 
with respect to u. The derivative with respect to u> leads to the vector Oi{Si-^.i) 
obeying the recursion relation 

with initial condition 

01(6^2) = J^^TiVi' exp(-/5ei +a;aiV?). (5.6) 
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At each step the vector Oi is rescaled as 



0^(5i+i) = [Oi(5,+i) - OiN',{Si^^)] /m 



{B.7) 



where rii and Oi are the sum of the elements of Ni{Si+i) and Oi{Si+i), respectively. 
The value oi q = q{ijj) is then obtained from the average of Oi 



In the limit of zero temperature one has to extract the diverging part in (4.8). 
This is achieved by means of a Legendre transform to express all the quantities in 
terms of the average energy per spin e = e(c<;, (3). The value of e for given 13 and u) can 
be obtained by introducing a vector of the derivative of Ni{SiJ^i) with respect to 
ad proceeding as done for q. 




1 



(S.8) 
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Figure captions 



Figure 1. Overlap probability P2{q;x) for x = 0.1 (a), x = 0.2 (b) and x = 0.5 (c). 
The dashed line is the average overlap probability P2{q) = Jq dxP2{q;x). 

Figure 2. Overlap probability distributions: P3{q; xi,X2) with impurities at xi = 1/3 
and X2 = 1/3 (full line), average probability Psiq) (dashed line) and gaussian 
distribution with mean value 1/3 and variance 2/27 (dot dashed line). 

Figure 3. Numerical results for the average overlap q as function of the temperature 
T = in the disordered system with p = 0.5 (a), p = 0.1 (b) and p = 0.01 (c). 

Figure 4. Numerical results for the average overlap q as function of the disorder 
concentration p at temperatures T = (a), T = 0.5 (b), T = 1 (c) and T = 2 (d). 

Figure 5. Spin glass susceptibility XSG = {{q ~ qY') ^ function of the temperature 
in the disordered system with p = 0.5 (a) and p = 0.1 (b). 
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